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ing. Also the influence of waveguide parameters to the second-harmonic 
spectra is addressed. 
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1. Introduction 

The process of second harmonic generation (SHG) has been one of the most frequently stud- 
ied nonlinear optical phenomena since the invention of lasers. This process was observed by 
Franken et al. in 1961 for the first time using crystalline quartz [1]. In principle, a second- 
harmonic (SH) field can be generated in many nonlinear materials that lack the inversion 
symmetry. However, phase matching of the nonlinearly interacting fields cannot be reached 
provided that the normal dispersion occurs [2] . Here, nonlinear anisotropic crystals with their 
birefringence have been found useful and allowed to fulfil phase-matching conditions. 

Alternatively, nonzero phase mismatch Ak among the interacting fields can be compen- 
sated by a periodic spatial modulation of the nonlinear susceptibility with period Apm (Apm = 
27r/|Ak|), as has been suggested by Armstrong 0. A commonly used technique in this case is 
based on inverting dielectric domains inside a ferro-electric nonlinear crystal using an intense 
static electric field flU . This process is called periodic poling (PP) and has allowed to utilize 
the highest material nonlinear coefficients (51 . 

The poling process leaves well defined domains of the alternating sings inside the nonlinear 
material. %^ nonlinearity can thus be approximated by a dichotomic function that can have 



significant coefficients at its first- as well as higher-order harmonics. Usually, a nonlinear struc- 
ture is designed such that only one harmonic is exploited. Contributions from other harmonics 
are considered as parasitic (71 in this case. If domain lengths are chosen such that several har- 
monics of nonlinear modulation are significant a rich spatial and spectral modal structure in a 
waveguide can be exploited O. This gives a substantial increase in the ability to tailor prop- 
erties of the nonlinear process inside a waveguide. Here, we study these abilities considering a 
waveguide that uses the first three harmonics of nonlinear modulation. 

Efficiency of the SHG process is inversely proportional to the transverse area of the fun- 
damental beam and, as a consequence, beam focusing may improve the efficiency of SHG. 
However, walk-off of the interacting fields puts limits to this property in bulk crystals. For this 
reason, waveguiding structures not suffering from this problem have become very important 
|9l . They have allowed to confine fields' energies into a very small area (typically of the order 
of tens of jlm 2 ) that results in the desirable high power densities |[T0l . Intense SH fields can 
then be obtained which, among others, allow to generate squeezed light (TTl[T2ll . 

Usually, waveguides support the propagation of the fundamental as well as several higher- 
order modes, for which intensity oscillations in the transverse plane are typical. The guided 
modes differ in their wave vectors (propagation constants). The phase-matching conditions and 
efficient SHG process can be reached simultaneously for several combinations of these spa- 
tial modes [Q31 [141 Q2) . For each of these individual phase-matched processes, certain spectral 
and also spatial characteristics of the fundamental beam are needed. We note that as we are 
interested in pulsed SHG we consider spectrally polychromanic modes determined by phase- 
matching conditions for spatially defined modes. Different spectral properties for different in- 
dividual processes are useful in tuning a specific nonlinear process. We note that, on the top 
of phase-matching conditions, also a sufficient spatial overlap of electric-field amplitudes of 
the interacting guided modes is needed. This disqualifies higher-order modes whose ampli- 
tudes change their signs several times in the transverse plane. Precise determination of efficient 
combinations of the modes is thus very important in characterizing the analyzed waveguide. 
This characterization is useful not only for the SHG process. Also the inverse process of spon- 
taneous parametric down-conversion well known for its ability to generate entangled photon 
pairs |[T4l [T5l [T6l benefits from the analysis. 

Waveguides are promising not only as efficient media for nonlinear interactions, they also 
allow for miniaturization and integration into more complex circuits. Their nonlinearities allow 
for the generation of squeezed light that is useful for optical communications and quantum 
metrology ifTTl . 

The theory of SHG in periodically poled waveguides is introduced in Sec. 2. Equations giving 
the transverse profiles of the fundamental and SH modes are derived in Sec. 3. The calculated 
and measured SH spectra are analyzed in Sec. 4. Spectral and spatial properties of type II SHG 
are discussed in Sec. 5. Conclusions are drawn in Sec. 6. 

2. Second-harmonic generation 

The evolution of electric-field amplitude E^ 2 ) of the SH field inside the waveguide is governed 
by the nonlinear wave equation (TTl US) • 



V = x^+y^ + zJ^. Linear permittivity tensor £ is assumed to have a diagonal form in the 
used cartesian coordinate system Constant c means speed of light in vacuum and /io 

is vacuum permeability. Vector P^ 2 ) of the second-order polarization introduced in Eq. ([]]) can 




(1) 



be expressed as p( 2 ) = 2£o d : E^E^, where d is the third-order tensor of nonlinear coeffi- 
cients and E^ denotes an electric-field amplitude of the fundamental pumping field. Symbol : 
stands for tensor shortening with respect to its indices. 

The electric field amplitudes of both the fundamental [E^] and SH [E^ 2 )] fields can be 
spectrally decomposed as follows: 



E^(x,y,z,t) = Y, L / d<o£$(z,(o)e$(x,y,G))exp (co)z- icot 



p=TE,TM n 



7 = 1,2. (2) 



The normalized electric-field mode functions e$(x,y,G)) form a basis and are obtained as 
eigenfunctions of the Helmholtz equation (see later in Sec. 0. The Helmholtz equation also 
provides the corresponding propagation constants f5pn(co). A detailed analysis of modal prop- 
erties of the studied waveguide has shown that two groups of modes differing in polarization 
exist |[T9l . Whereas (quasi-) TE modes have electric-field polarization vectors nearly parallel to 
the x axis, electric-field polarization vectors of (quasi-) TM modes are nearly perpendicular to 
the x axis. In addition, there exists a system of transverse modes labelled by index n. Finally, 
envelopes &pn(z, co) of the spectral amplitudes of the corresponding modes propagating along 
the -f z axis have been introduced in Eq. ©. 

Substituting the decomposition of electric-field amplitudes E ^ in Eq. © into the wave equa- 
tion (OQ) a set of nonlinear differential equations for the envelope amplitudes Spi} (z, co) can be 
derived. Invoking the slowly varying amplitude approximation and non-depleted fundamental 

(2) 

pump field approximation, the following formula for envelope amplitude at the end of the 
waveguide of length L has been obtained after simple integration: 

2 

4 2) (L,ft>2) = tft>2 £ ZfdmDfcfa^ltttL, a*, on) 

J) q(2)/^ \ fc c =TE,TM Lm J 



x4%tti)4^(0,a)2-coi). (3) 

We note that the integration over frequency (0\ incorporates all frequency contributions occur- 
ring in the polychromatic fundamental field. In Eq. ©, effective nonlinear coefficients 
have been introduced: 

D a k f^(co 2 ,co l ) = d M I dx [ JyVie^fott (4) 

J OO J OO 

The effective nonlinear coefficients Dj^J£ incorporate in their definitions both polarization prop- 
erties and overlap of the transverse parts of the interacting modes amplitudes. Constant dM gives 
the amplitude of decomposition of the actual modulation of nonlinear coefficient into the M-th 
harmonic. The coefficients Dff^ are assumed to be weakly frequency dependent. Considering 
nonlinear coefficients Dff^ involving higher-order spatial modes, their values are usually small 
due to frequent changes in the sign of the electric-field amplitudes of these modes in the trans- 
verse plane. The coupling functions r^(L, COz, CO\) introduced in Eq. © characterize fields' 
evolution along the z axis: 

The nonlinear phase mismatch AjS^ occurring in Eq. (0) is defined as 

APt ®i) = P%\o>i) - p£ ] (o>i) - /3« (0)2 - 0)x) - (6) 




Fig. 1. Topo graphs of (a) experimental and (b) theoretical intensity / = | 2 in the 
transverse plane formed by white light propagating inside the waveguide. Also cuts along 
the lines indicated in the topo graphs are shown. 



the last contribution 2;tM/Apm originates in the M-th harmonic of the periodic modulation of 
nonlinear coefficient. 

3. Equations giving fields' spatial profiles 

Material of the waveguide, KTP, is anisotropic and is positioned such that only diagonal ele- 
ments £ x , e y , and £ z of the linear susceptibility tensor are nonzero. The crystallographic axes 
of KTP (x c , y c , z c ) coincide with the axes of the coordinate system (z, x, y). The elements £ x , 
Ey, and £ z can be expressed using indices of refraction (£^ = t; = x,y,z, that take the 
following form in the studied waveguide l9l l20ll2H : 

H \ y >o = n ^ + An ^ ct [-w/2,w/2] 0) erfc (-y/h) , 

n s( x >y)\ y <o = 1; Z=x,y,z- ( 7 ) 

The indices of refraction n^ characterize the substrate. Changes An% of the indices of refrac- 
tion are introduced in the process of waveguide formation. For X\ = 800 nm (A-2 = 400 nm) 
the following values have been used: n x = 1.75719, An x = 0.009, n y = 1.84546, An y = 0.013 
(n x = 1.84435, An x = 0.018, n y = 1.96775, An y = 0.019) (20|. Above the waveguide, air is 
assumed. Waveguide's depth is characterized by parameter h, whereas parameter w denotes 
waveguide's width. We have assumed h = 10 /im and w = 5 /im for the analyzed waveguide. 
These values have been checked experimentally by measuring the profile of white light leav- 
ing the waveguide and comparing the obtained profile with that arising from the model (see 
Fig.[T]). The rectangular function rect[_ w / 2 ,w/2] (x) equals 1 in the interval (— w/2,w/2) and is 
zero otherwise. Also the error function erfc(x) has been invoked in Eq. ©. 

Usually, TE modes are conveniently obtained by solving the wave equation for the electric- 
field amplitude. The wave equation for the magnetic-field amplitude is then useful in deter- 
mining TM modes. However, the spatial dependence of indices of refraction in the studied 
waveguide causes difficulties and asymmetry in the obtained equations. As a consequence, the 
use of wave equation for the magnetic-field amplitude H, 

Vx(^VxH) = -'^, (8) 



is preferred in both cases. 

The substitution H(x,y,z,t) = h(jc, y, (o) exp(/j3z — i(Ot) in Eq. ([8]) provides a set of equations 
for cartesian components of the magnetic-field envelope h(x,y, co) of a mode characterized by 
frequency CO and propagation constant j3 . Considering the relation 

h z = i(dh x /dx + dh y /dy)/p (9) 

originating in the Maxwell equation VH = 0, two coupled differential equations for the com- 
ponents h x and can be derived: 



-*o)h*, (10) 
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ko = (o/c. The last component h z (x,y,a>) is then obtained from Eq. (0 provided that the 
remaining components have been determined. Knowing the magnetic-field envelope h, the 
electric-field envelope e is derived from the Maxwell equation that gives e(x,y, co) = V" -1 [/V x 
h(x,y, co) — jSz x h(x,y, o))]/(£o<w); z means the unit vector along the z axis. 
Considering Eqs. (fTOt and (fTTT) . the cross contributions 
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are usually very small and so they can be omitted [22] . As a result, equations (ITT)b and (fTTT) 
become independent and so can be solved separately giving (quasi-) TE and TM modes. This 
considerably simplifies the analysis. The obtained equations can be written as follows: 



d 2 h x 

dx 2 + 
£xd 2 h y | 
£ z dx 2 dy 2 



£ y d 2 h x 
£ z dy 2 
d 2 h v 



(l3 2 -£yk 2 )h x , 
(l3 2 -£ X k 2 )hy. 



(12) 



For certain profiles of permittivity £ , the solution of Eqs. (1121) can be found analytically 
E21 (231 . For example, the simplest approximation of waveguide's profile by the rectangular 
functions in both x and y directions together with the known analytical form of eigenmodes has 
been used in lfT4ll to interpret the experimental results. However, only the numerical approach 
is possible for transverse profiles of real waveguides. For this reason, a large variety of sophisti- 
cated numerical methods has been developed l22l . Among them, a finite difference method has 
become popular because of simple implementation, even in its full- vector form. Unfortunately, 
it requires a large amount of computer resources for solving the eigenvalue problem. This rep- 
resents a problem when more complex transverse profiles are considered. Nevertheless, it has 
been successfully applied both in its scalar l24ll as well as semi- vectorial (l3j[25l forms. 

On the other hand, a finite elements method (22l [261 has been found more suitable due to 
its stability and lower demands on computer resources. It is resistant against non-physical so- 
lutions. Moreover, the use of higher-order approximation functions results in lowering of the 
dimension of eigenvalue problem. For this reason, many commercial waveguide mode-solvers 
are based upon this method. For instance, it has been applied for the determination of waveg- 
uide's mode dispersion in (27). In our investigation, we have developed an implementation of 
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Fig. 2. Topo graphs of calculated electric-field intensity distributions 7 tr (7 tr (x,y) = 
|e(x,y) | 2 ) in the transverse plane for the first three spatial modes denoted as (0,0), (0,1), 
and (1,0) for the fundamental TE-polarized (a) and SH TE-polarized (b) field; w = 5 /im, 
h = 10 /im. Frame boxes in the graphs indicate boundaries of the waveguide. 

the finite elements method that is based on the Galerkin method ll22l . Easy parallelization of 
this implementation makes it superior above other approaches. 

Intensity profiles of the first three modes both for the fundamental and SH TE-polarized fields 
as they emerge from the developed method are depicted in Fig. [2] The transverse profiles of 
TM-polarized modes are similar to those appropriate to the corresponding TE-polarized modes. 
According to the convention, spatial modes are labelled by two numbers that give the number of 
nodes along the x and y axes. The numerical analysis has revealed that the analyzed waveguide 
supports more than 30 spatial modes at X2 = 400 nm for both TE and TM polarizations and 3 
(5) spatial modes at X\ = 800 nm for TE (TM) polarization. 

4. Experimental second harmonic generation in the waveguide 

The analyzed waveguide was fabricated on a KTP substrate (10. 5x2x1 mm 3 ) by Rb + ion 
diffusion together with other cca 50 similar waveguides (manufacturer AdvR Inc.). The re- 
quirement for an efficient type II SHG at X\ = 800 nm has resulted in poling the nonlinear 
material with period Apm = 7.62 /im. In the fabrication process both horizontal and vertical 
dimensions of the waveguide are under control. The horizontal width is given by geometry 
of the lithographic mask that leads to a profile with well defined boundaries owing to highly 
anisotropic diffusion. On the other hand, the profile in the vertical direction is formed by ion 
diffusion gradually in time. At some instant the profile of refractive index becomes temporally 
stable and can be approximated by the error function lETTl . 

Due to imperfections in the fabrication process, fluctuations in waveguide's parameters in- 
cluding its width w, depth h, duty ratio of the poling period as well as the poling period itself 
inevitably occur. Moreover, there exist several areas without poling inside the structure and 
also the poling duty ratio may change from 50 % to 75 %. These fluctuations result in broad- 
ening of the SH spectrum. To judge the amount of spectral broadening we have plotted the 
dependence of SH wavelength shift AX2 on waveguide's width w, depth h and poling period 
Apm in Fig. [3] for five individual nonlinear processes dominating in SH spectrum. According 
to the manufacturer [28|, the expected variations in waveguide's width w, depth h and poling 
period Apm are in turn ±0.1 /im, ±2 /im and ±0.1 /im. This in accord with the curves of 
Fig. [3] results in the estimated spectral SH widths equal approximately to 0.05 nm, 0.1 nm, and 
2 nm, respectively. Fluctuations in the poling period Apm thus represent the main source of SH 
spectral broadening observed with a sufficiently broad fundamental spectrum. Imperfections 
in nonlinear modulation even lead to the occurrence of three SH fields originating in different 
polarization configurations and having different ideal values of poling period Apm (see below). 



Larger differences in the values of width w, depth h and poling period Apm lead to significant 
shifts of spectral positions of individual lines that result in the change of the overall SH spectral 
profile. The character of this change is illustrated in Fig. H showing the calculated SH intensity 
spectrum / for three values of the width w. 




Fig. 3. SH wavelength shift AA2 as it depends on waveguide's width w (a), depth h (b) and 
poling period Apm (c) for the following individual processes: (0,0) + (0,0) — > (0,0) (black 
curve), (0,0) + (0,0) (0,1) (red curve), (0,1) + (0,1) (0,0) (green curve), (0,1) + (0,1) 
— >> (0,1) (blue curve), and (1,0) + (0,0) — > (1,0) (magenta curve). In (c), all curves nearly 
coincide. 



To observe the full profile of a spectral line characterized by its width AA2, the width AAi 
of the fundamental field has to be sufficiently large. Consideration of energy conservation in 
the nonlinear process provides in this case the relation AAi = 2V2AA2. However, the needed 
width AAi has to be larger than 2\flbJi2 because of phase-matching conditions that put addi- 
tional constraints to the nonlinear process. This relation between the widths AAi and AA2 also 
determines how wide area in the SH spectrum is observed for a given fundamental spectral 
width AAi • We note that a more detailed information about the phase-matching conditions and 
their spectral dependence inside a broadened line can be reached by sum-frequency generation 
utilizing spectrally narrow fields fT3l . 

The experimental investigation of SHG process has been carried out in the setup sketched 
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Fig. 4. SH spectral intensity / depending on wavelength A2 and calculated for different 
values of waveguide's width w; 1(002) = Hak \&ak{L, &>2)| 2 - Solid curves correspond to 
ideal spectral resolution 0.13 nm (given by phase-matching conditions along the z axis 
for a 10.5 mm long waveguide with the inclusion of material and waveguiding dispersion). 
Dashed curves are appropriate for the expected real spectral broadening of 1 nm; w = 4 jim 
(blue curve), w = 5 /im (red curve) and w = 6 /im (green curve), h = 10 jUm. 
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Fig. 5. Sketch of the experimental setup for waveguide (WG) investigation. Symbol Ft 
denotes interference filter, Fa neutral density filter, Fg edge-pass filter, A/2 half- wave plate, 
LP linear polarizer, O microscope objective, DB dichroic beam-splitter, and CCD means 
CCD camera. 



in Fig. [5] Light from a tunable femtosecond Ti-sapphire laser (87 MHz repetition rate, 100 fs 
pulse duration, central wavelengths in the range 790 — 810 nm, spectral width 10 nm) has been 
attenuated combining half- wave plate (A/2) and linear polarizer to attain the power of about 
20 mW. Another half- wave plate (A /2) has been used to control the pump-beam polarization. 
The pump-beam central wavelength has been adjusted by tilting an interference filter defining 
the spectral width of 3 nm. A small fraction of the beam has been deflected to allow for precise 
monitoring of the pump-beam position in the transverse plane needed for efficient coupling of 
the beam into the waveguide using a 10 x microscope objective. The light leaving the waveguide 
and composed of both the pump and SH beams has been imaged by a 20 x objective with a 
larger numerical aperture. A dichroic beam-splitter has separated the pump beam from the SH 
beam. The obtained beams have been analyzed by power-meter, CCD camera or spectrometer 
(resolution 0.1 nm) after passing through a linear polarizer. As the indices of refraction in KTP 
are temperature dependent, the KTP chip has been thermally stabilized at 25 °C using a Peltier 
element with temperature sensor and feedback loop. We note that the ideal phase-matching 
condition shifts by cca 1 nm per 25 °C. 

For the chosen orientation of KTP, three types of SHG processes differing in fields' polariza- 
tion properties can occur simultaneously due to different harmonics of the periodical modula- 
tion of nonlinearity. In type SHG, both the fundamental and SH fields have TM polarizations. 
If the fundamental beam is TE polarized and the SH beam TM polarized type I SHG is ob- 
served. A TE-polarized SH beam occurs only if both TE- and TM-polarized components of 
the fundamental beam are present (type II SHG). Different conversion efficiencies are found 
for the discussed processes because of different values of used nonlinear coefficients, overlaps 
of mode functions, and also phase-matching conditions. Especially, phase-matching conditions 
depending on the value of poling period Apm allow to control these efficiencies. The analyzed 
waveguide was primarily optimized for type II SHG giving Apm = 7.62 jim. Nevertheless, 
type SHG (with the optimum poling period 3.08 /im) can be observed exploiting the second 
harmonic of nonlinear periodic modulation. The presence of second harmonic of nonlinear pe- 
riodic modulation together with all even harmonics is caused by the declination of the actual 
nonlinear modulation from the ideal shape with positive and negative domains of equal lengths 
l29l . We note that a nonlinear modulation with equal lengths of positive and negative domains 
has only odd harmonics nonzero. Also, the third harmonic of nonlinear periodic modulation 
allows to approximately arrive at the phase-matching conditions of type I SHG for which the 
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Fig. 6. Measured normalized power P of SH field as it depends on angle a giving the ori- 
entation of linear polarization of the incident pump beam for (a) TE- and (b) TM-polarized 
SH field (a = deg corresponds to TE polarization). Curves characterizing three different 
waveguides are shown for comparison to judge fabrication capabilities. The curves were 
obtained with experimental relative error 4 %. 



poling period 1.83 /im would have been optimum. The experimental nonlinear conversion ef- 
ficiencies r/ in this configuration are summarized in Tab. [T] Type II SHG is the most efficient 
process due to optimum phase-matching conditions and despite the lower value of nonlinear co- 
efficient J32 compared to that of coefficient used in type SHG. Also better phase-matching 
conditions of type I SHG allow more efficient SHG than in type 0. However, and importantly, all 
three types of SHG have comparable efficiencies. Their relative values can suitably be changed 
by changing the poling period Apm- 





Process 


d [pm-V- 1 ] 


77 [W^-cnr 2 ] 


type 


TM + TM -> TM 


d 3 s = 10.7 


2.83 ± 0.06 


type I 


TE + TE ^ TM 


J32 = 2.65 


4.44 ±0.12 


type II 


TE + TM ^ TE 


J32 = 2.65 


4.78 ± 0.03 



Table 1. Types of observed SHG processes, their polarizations, used nonlinear coefficients d 
(T8), and experimental nonlinear conversion efficiencies 77 are given. Conversion efficiency 
77 has been determined as r\ = PsugP\^ 2 L~ 2 , where Pshg C^in) gives the power of outgoing 
SH (coupled incident pump) beam and L is the length of nonlinear medium ifTQl . 

Suitable choice of polarization directions of both the fundamental and SH fields allows to 
extract individual types of SHG. The intensity of SH field in the most efficient type II SHG 
process can be separated from the overall field by using a polarization analyzer aligned along 
TE polarization. The experimental evidence of observing this process is given by measuring 
periodic oscillations of SH intensity depending on the orientation of linear polarization of the 
incident pump beam [see Fig.Oa)]. Aligning the polarization analyzer in SH field along TM 
polarization the remaining two processes are visible. Changing the orientation of linear polar- 
ization of the incident pump beam, we continuously move from one process to the other [see 
Fig. Ob)]. Whereas TM polarization of the pump beam supports type SHG, TE polarization 
is suitable for type I SHG. 

We further pay attention to the most efficient type II SHG process. The necessity to have 
both TE- and TM-polarized fundamental field distinguishes it from the remaining two types. 
This is important from the point of view of the reverse process of parametric down-conversion 
[1741 1731 1771 l30l l3T1l in which pairs of photons with orthogonal polarizations are emitted. The 
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Fig. 7. Experimental SH spectral intensity / depending on wavelength A2 and topo graphs of 
spatial electric-field intensity profiles 7 tr (x, y) of TE-polarized pump beam, TM-polarized 
pump beam, and TE-polarized SH beam are shown in sequence for individual processes: 
(a) (0,0) + (0,0) -> (0,0), (b) (0,1) + (0,1) -> (0,1), and (c) (1,0) + (0,1) -> (1,0). 



orthogonal polarizations of two photons would then allow to spatially separate both photons 
using a polarizing beam-splitter. 

5. Spatial and spectral properties of type II second harmonic generation 

In agreement with the numerical analysis, the experiment has revealed three (five) spatial TE- 
(TM-) polarized modes [(0,0), (0,1), (1,0), (1,1), and (0,2)] in the fundamental field inside the 
waveguide (see Fig. [7]). As expected, the spatial modal structure of the SH field has been found 
richer owing to doubled relative waveguide dimensions with respect to the SH wavelength. 
The experimental spatial intensity profiles shown in Fig. [7] are partially distorted in comparison 
with the theoretical ones. This is caused by waveguide imperfections, non-ideal coupling of the 
fundamental beam leading to the excitation of more than one spatial mode and subsequent SHG 
into several spatial modes. 

In the nonlinear process, a sufficient overlap of three chosen spatial mode profiles is needed to 
observe an individual nonlinear process in the measured SH spectrum. As also phase-matching 
conditions play an important role, a given individual process occurs only in a restricted spectral 
range of SH field (see Fig. [7]). If spatial modes of the fundamental field are uniformly excited 
by a spectrally broad pumping, several individual nonlinear processes typically occur in the 
waveguide and contribute to the SH spectrum. As a consequence, the expected SH intensity 
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Fig. 8. Theoretical SH spectral intensity / observed in an ideal (red curve) and real (blue 
curve) structure. Individual combinations of spatial modes are identified. 

spectrum is relatively broad (see Fig. [8]). As SH spectra belonging to individual processes are 
relatively wide and close to each other, only spectral filtering of the SH field itself is not suffi- 
cient to distinguish them. After spectral filtering, individual processes become experimentally 
available by exciting only specific spatial modes in the fundamental beam. Technically, the 
proper excitation of the chosen spatial mode(s) is accomplished by suitable coupling of the 
pump beam into the waveguide. 



Modes 




^2,exp 


^2,theor 


A^2,exp 


^exp 


Aheor 


(TE + TM - 


► TE) 


[nm] 


[nm] 


[nm] 


[arb.u.] 


[arb.u.] 


00 + 00 ^ 


00 


398.3 


398.6 


0.8 


1.00 


1.00 


00 + 00 -+ 


01 


394.2 


394.4 


0.9 


0.65 


0.81 


10 + 00 ^ 


10 


398.8 


399.5 


0.6 


0.42 


0.78 


01 +01 ^ 


00 


402.8 


403.0 


0.8 


0.38 


0.43 


10 + 01 ^ 


10 


401.1 


401.9 


0.9 


0.35 


0.42 


01 +01 -+ 


01 


398.7 


397.1 


0.6 


0.32 


0.18 



Table 2. Individual SHG processes identified in the waveguide, experimental central SH 
wavelength A2, e xp> theoretical central SH wavelength A2 ? theor> spectral width AA2, e xp> exper- 
imental relative power Pexp? and theoretical relative power ^theor ^e given. Experimental 
central wavelengths are measured with the error 0. 1 nm given by the spectrometer response 
function. Spectral widths (relative powers) have the relative error 5% (9%). 

Careful alignment of the pump-beam coupling and use of 3 nm wide spectral filters in this 
beam have allowed to observe cca 10 spectral lines from the SH spectrum plotted in Fig. [51 SH 
central wavelengths A2, spectral widths AA2 as well as relative powers of the observed SHG 
processes are summarized in Table [2] The comparison of theoretical and experimental param- 
eters in Table [2] has revealed good agreement in central wavelengths of the SH lines. Larger 
differences between the theoretical and experimental relative powers P can be explained by the 
difficulties in reaching an efficient coupling of the corresponding pump mode. The experimen- 
tal spatial intensity profiles of the three interacting modes and the corresponding SH intensity 
spectra are plotted in Fig. [7] for three typical mode combinations. These graphs document real 
capabilities in controlling spatial properties of the SH beam. If the analyzed waveguide were 



used for parametric down-conversion, several individual processes based on different spatial 
modes would simultaneously occur contributing together to the signal and idler fields' profiles 
[e.g., parametric processes (0,0) -> (0,0) + (0,0) and (0,0) (0,1) + (0,1)]. This might be used 
for the generation of modally entangled photon pairs ll32l . We finally note that temperature 
tuning is important in the determination of efficient individual processes. This gives additional 
possibilities in controlling the properties of the generated fields. 

6. Conclusions 

The process of second harmonic generation in a periodically-poled KTP waveguide has been in- 
vestigated both experimentally and theoretically. Three types of the nonlinear processes (type II, 
I, and 0) have been observed simultaneously utilizing the first, second and third harmonics of 
the spatial nonlinear modulation. A scalar finite elements method based on the Galerkin method 
has been applied to calculate spatial mode profiles, propagation constants, and frequencies of 
the interacting fields. It has been shown that positions of the spectral lines of individual pro- 
cesses depend strongly on waveguide's depth and width in ideal waveguides. In real waveg- 
uides, fabrication imperfections considerably broaden single spectral lines occurring in ideal 
waveguides. Individual spatial and spectral processes in type II interaction have been experi- 
mentally characterized by changing pump-beam spectral filtering and coupling into the waveg- 
uide. The analyzed waveguide has been recognized as a versatile source of second-harmonic 
light. 
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